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We present a bond-operator mean field theory for the 
Kondo lattice model at half filling in two (2D) and three 
(3D) dimensions. A continuous quantum phase transition 
from an antiferromagnetic to a spin-gapped singlet ground 
state is found at J/t = 1.505 (1.833) in 2D {3D). Addition- 
ally we evaluate the quasiparticle dispersions as well as the 
staggered magnetic moment and provide a comparison with 
complementary numerical approaches. 
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The Kondo Lattice Model (KLM) describes the 
exchange-scattering of a band of itinerant conduction- 
electrons at a lattice of localized magnetic moments. It 
serves as a basic model for heavy fermion materials in 
the integral- valent limit ^ . At half filling of the conduc- 
tion band it is believed to give a description of Kondo 
insulators, which have been of considerable interest in 
recent years [Q. In one dimension the ground state at 
half filling is a spin singlet for all values of exchange- 
scattering and conduction-band widths Q . In higher di- 
mensions it has been suggested early on, that the com- 
petition between Kondo screening and the Ruderman- 
Kittel-Kasuya-Yoshida (RKKY) interaction leads to a 
quantum phase transition between a global spin-singlet 
and an antiferromagnetically ordered phase |^,^. This 
scenario has been corroborated in two dimensions by vari- 
ational 1^ calculations, series expansion [Q and mean- 
field H] approaches. Numerically exact results in 2D have 
been obtained recently by QMC . In 3D only series 
expansion is available Q . 

The purpose of this work is to introduce a novel mean- 
field theory for the KLM at half filling in two and three 
dimensions. In contrast to other mean-field calculations 
1^,^ our approach is based on a bond-operator represen- 
tation of the KLM which is suitable for strong exchange 
scattering and has proven to be useful in dimerized spin 
systems p2| . Moreover, our treatment goes beyond re- 
cent mean-field work focusing on the Kondo-necklace 
problem which neglects conduction-electron charge fluc- 
tuations [HI. The KLM reads 



H = -t , 



(1) 



with spin operators Sc(/) = \ Y.a,a' c{f)lT„,a'c{f)^, and 
destruction(creation) operators cl'l and f^^} for itinerant 



c and localized / electrons of spin a at sites i. An f- 
occupation of exactly one per lattice site is implied. 



The local Hilbert space consists of one / electron with 
spin up or down and additionally up to two itinerant elec- 
trons. The resulting eight possible states can be created 
by applying the following operators onto the vacuum |0) 
of an empty site 
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where the operators s and t are equivalent to the so- 
called bond operators of and are assumed to obey 
bosonic commutation relations. The fermionic operators 
a and b have been introduced first in and label states 
with one or three electrons per site. In order to suppress 
unphysical states a constraint of no double occupancy 

+ E ^Jo., + E + E ^.A, = 1 (3) 



has to be fulfilled. The original fermion and spin opera- 
tors are represented by 
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with Sa(b),j = kT>a,a''^ib)l,j'r'y,<y'a{b)^,^j and p| = 1, 
Pl = — 1. This mapping is exact and preserves the proper 
commutation relations provided that s and t are bosons, 
a and b are fermions, and (|^) is satisfied. 

Rewriting the KLM in terms of (||,§-^) leads to 
a strongly correlated boson-fermion model which can- 
not be solved exactly, but approximations have to be 
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made. Analogous models in the context of quantum- 
disordered ground states in dimerized spin systems have 
recently been treated by various schemes of approxima- 
tion |l4|-|l|]. Here we are interested in a mean-field 
description of the transition from an antiferromagnetic 
state to a spin-singlet regime. The latter can be described 
by allowing for a condensate of singlets 



(7) 



while the antiferromagnetically ordered phase requires an 
additional condensation of one of the triplets 



(tz.,) = (4,,) = = i-iym 



(8) 



For the remainder of this work we assume a bipartite lat- 
tice structure with the factor of (—I)-' being a shorthand 
for '-hl(-l)' on the A(B) sublattice. 

Inserting (|^,^) into and (|l|) we drop all longitudi- 
nal and transverse spin-fluctuations, i.e., we only keep the 
mean-field values of s and t in addition to the fermions 
a and b. We obtain 
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The first sum cx t describes hopping of particle- or hole- 
like fluctuations with a dispersion renormalized by a fac- 
tor of (s^ — m?), taking into account the alternation 
of the sign of m^. This prefactor leads to a reduction 
of the fermionic kinetic energy in the antiferromagnetic 
phase. The second sum oc t generates the mean- field ana- 
log of the intermediate state of the RKKY process and 
its hermitian conjugate, i.e., the destruction(creation) of 
two adjacent two-particle states accompanied by the cre- 
ation(destruction) of a pair of two adjacent three- and 
one-electron states. Second order processes from this 
term drive the transition into the magnetic state. Thus 
we believe that this mean-field Hamiltonian incorporates 
the basic ingredients to induce the quantum-phase transi- 
tion of the KLM. Furthermore in (||) we have introduced 
the usual chemical potential A to set the global particle 
density and a local Lagrange multiplier /i^ in order to to 
enforce the constraint M). 




FIG. 1. Mean-field order parameter s, m and the staggered 
moments Mc and Alf as function of J/t for the 2D KLM. 



To diagonalize the Hamiltonian we first replace the lo- 
cal Lagrange multiplier by a global one fii /i. Sec- 
ond we switch to an antiferromagnetic unit cell and fi- 
nally perform a generalized Bogoliubov transformation 
to eliminate particle- number non-conserving term of type 
ab{a^b^). This leads to 4 bands wi,2(k) = X±Ei{k) and 
a;3,4(k) = A ± £'2(k) which are twofold degenerate by 
spin-z quantum number 



E2M = \//^' + oCkC™" + + 2VFk 
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Here = ^'^tJ2d=i cosfc^. At half filling the lower two 
bands, i.e. tj^2,4, are completely filled while the upper two 
bands , i.e. ti^i,3, are empty. This leads to a ground state 
energy of 
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which is independent of A. 

The staggered magnetization of the c(/) elec- 

trons is obtained from a direct evaluation of the corre- 
sponding matrix elements in the mean-field ground state 
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The second term in Mf is due to a staggering of the 
spin density of the a and b fermions which develops at 
non-zero values of m due to the the second sum oc i in 
(^. This terms induces a spin-dependent non-diagonal 
momentum-space component in the a and b Greens func- 
tions at the antiferromagnetic nesting wave- vector. 
The mean-field equations 
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have to be solved self-consistently for the order param- 
eters s, m and In the magnetic phase (m ^ 0) the 
system of equations (pl5[) can be written as 



= 2J 



1 

2N ^ 

k 



E2M Eiu 



= -J + 4At--^5]2e2(m2 + ,2) 



2N 

2N ^ 



k 

2 I o2^3 



= s^ 



1 - 



2W^k 

-Em 



£^2,k 
2 



2A^ 



-B2,k 

2 

£^2,k 
2 



2 



2iV ^ 

k 



4W^k 



2,k 



2 



(14) 



while for the disordered one (to = 0) it simplifies consid- 
erably 
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For Ek = 0, i.e. t = 0, only (|T^) has a solution. This 
solution also provides for the correct quasiparticle gap, 
i.e. fi = jJ. Upon increasing t charge fluctuations, 
i.e. creation of a- and 6-excitations contribute to the 
ground state reducing the stability of the singlet state. 
At J/t the ground state stems from the solution of 
( p^ ) with s = m. This state displays a complete polariza- 
tion of the / spins, while the c spin density is polarized 
only partially. The latter is consistent with the itinerant 
character of the c fermions, implying a finite density of 
empty and doubly occupied sites, i.e. a finite density of a 
and b fermions. Since s = to at J/t = the diagonal part 
of the kinetic energy of the a and b particles vanishes at 
this point. 

At intermediate J/t we determine the ground state by 
solving ( p^ , p^ ) numerically. Fig. ^ shows the singlet and 
triplet order parameters as well as the staggered magneti- 
zations for the 2D KLM. We observe a continuous quan- 
tum phase transition from the singlet to the magnetic 
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FIG. 2. Quasiparticle dispersions for (a) J/t=1.2 and (b) 
J/t=2. 



phase at {J/t)c — 1.505. This is in good agreement with 
the data of a recent QMC study |^ which has determined 
the phase transition to occur at {J/t)c — 1.45±0.05. Sim- 
ilar values of {J/t)c have also been reported from vari- 
ational Monte-Carlo simulations (|] ((J/t)c = 1.4 ±0.1) 
and series expansion [|((J/t)^ = 1.43 ± 0.2). From Fig. 
^ it can be seen that the maximum magnetization, i.e. 
s = to, prevails only at J = with a continuous increase 
of s vs. TO, i.e. screening of the local moment, to occur 
as J approaches Jc- Such coexistence of Kondo screen- 
ing and antiferromagnetic order has also been reported 
recently from QMC calculations [0 for all J < Jc and 
within a small window of values of J from a mean-field 
study B. 

Fig. l2 shows the quasiparticle dispersion of the occu- 
pied bands for two values of J/t which are in the sin- 
glet and the magnetic phase, i.e. (a) J/t = 2 and (b) 
J/t = 1.2 respectively. These bands are split by a gap 
from the unoccupied bands which are located symmetri- 
cally reflected along the line E/t = at positive energies. 
For J > Jc the four bands a;i...4(k) collapse onto only two 
bands by a mere backfolding which has been carried out 
in fig. ^(a) leaving a single occupied band to be displayed. 
For J < Jc two distinct bands are present throughout the 
entire magnetic Brillouin zone. Since the Hamiltonian 
(^ incorporates scattering with a magnetic wave-vector 
k = (tt, tt) in the off-diagonal terms only, no additional 
gap opens along the line k^r with fcjv.a; + fcjv,^ — ""i 
W^kw = in (Eoh. In this context we note, that the 
interpretation of the band-gap in this mean-field theory 
changes quasi- continuously from a gap induced by singlet 
formation at J > Jc to a magnetic gap as J — > 0. 

In fig. ^ we display a set of results identical to that of 
fig. |] however for the 3D case. Here the phase transition 
occurs at a slightly larger value of {J/t)c = 1.833, which 
is in reasonable agreement with an estimate of {J /t)c — 
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FIG. 3. Mean-field order parameter s, m and the staggered 
moments Mc and Mf as function of J/t for the 3-D KLM. 

2.04±0.16 reported from series expansion Q. Again, the 
/ spins are fully polarized as J — > 0, while the maximum 
value of Mc is nearly identical to that of the 2D case. 
To conclude several comments are in order. First, and 
much in contrast to usual approaches to the KLM 
our method is best suited for the limit of strong 
and quasi local Kondo screening at large J/t. In that 
limit the Kondo effect can be viewed as a molecular sin- 
glet formation within each unit cell resulting in an al- 
gebraic energy scale of order J, rather than the usual 
Kondo energy-scale Tk ~ iexp(— i/J). While the large- 
J limit may obliterate some of the subtleties genuine to 
the Kondo effect at J ^ t, we believe that it is a superior 
starting point for analytic studies of the quantum phase 
transition in the 2D and 3D KLM since this transition 
occurs at J/t > 1. Second, we note that while we have 
neglected quantum fluctuations of the triplet order pa- 
rameter, it would be interesting to incorporate them into 
future studies. In particular, it is conceivable that trans- 
verse fluctuations due to the t^^^y-operators will reduce 
the staggered magnetization in the ordered phase. This 
seems consistent with QMC finding a smaller magnetiza- 
tion 1^ than we observe within the mean-field approach. 
Finally an extension of the scheme presented here to in- 
corporate Coulomb correlations or finite doping, off half 
filling, into the conduction band are open issues. 

In summary we have studied the KLM using a novel 
bond-operator mean-field theory. In good agreement 
with complementary approaches we find a quantum 
phase transition at {J/t)c = 1.505(1.833) in 2(3) dimen- 
sions. In addition we have evaluated the magnetization 
in the ordered phase and the quasiparticle dispersions. 
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